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1. INTRODUCTION 


The task of this research, as originally proposed, was to focus on 
the development of a statistical methodology that may be used to enhance 
the Space Shuttle Main Engine (SSME) performance prediction accuracy. 
This methodology was expected to be used in conjunction with existing 
SSME performance prediction computer codes to improve parameter 
prediction accuracy and to quantify that accuracy. However, after a 
review of related literature, the researchers for this project concluded 
that the proposed problem required a coverage of areas such as linear 
and nonlinear system theory, measurement theory, statistics, and 
stochastic estimation . 

Since state space theory is the foundation for a more complete study 
of each of the before mentioned areas, these researchers chose to 
refocus the emphasis of the project to cover the more specialized topic of 
state vector estimation procedures. State vector estimation was also 
selected because of current and future concerns by NASA for SSME 
performance evaluation. That is, there is a current interest in an 

improved evaluation procedure for actual SSME post flight performance as 
well as for post static test performance of a single SSME. Furthermore, 
Taniguchi (1985) and Cikanek (1985) described a current investigation of 
analytical methods that may be used to improve test stand failure 
detection. Many of the analytical tools that have been suggested for 
failure detection applications are based on a function of a state vector 
estimate. Hence, this paper considers the issue of post flight /test state 
variable reconstruction through the application of observations made on 
the output of the Space Shuttle propulsion system. The incentive for 
the concepts discussed here is rooted in a NASA contract report by 
Rogers (1987). Rogers used the 
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Kalman filtering procedure to reconstruct the state variables of the Space 
Shuttle propulsion system. 

An objective of this paper is to give the general setup of the Kalman 
filter and its connection to linear regression. A second objective is to 
examine the Rogers (1987) reconstruction methodology for application 
to the reconstruction of the state vector of a single Space Shuttle Main 
Engine (SSME) by using static test firing data. 

Throughout this paper, underlined capital letters are used to denote 
vectors, capital letters denote matrices, and the identity matrix is denoted 
by the letter I. All vectors are of the column type and the transpose 
of any matrix or vector is denoted by using the letter T at the superscript 
position. The letter E denotes the expectation operator, N(U,E) denotes 
the multivariate normal probability distribution with mean vector U and 
covariance matrix I , The caret symbol " ~ ” written directly above -a scalar 
or vector denotes a statistical estimator of that vector or scalar. 

2. ESTABLISHMENT OF NOTATION 

For any Space Shuttle flight, let Z t denote the observed values of 
vector Z at time t. Each component of Z t represents a relevant ouptut 
of the Space Shuttle propulsion system that can be measured. For example, 
the components of Z t may be chamber pressure, oxygen flow rate, hydrogen 
flow rate, etc., for each of the three SSME’s. Rogers (1987) lists 35 
components for vector Z t and 71 components for vector X t where X t is 
defined below. 

Vector X t is a state vector of parameters to be estimated at time 
t. It is assumed that the observation vector Z is a function of the 
state vector X t . That is, 
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where h is some function and -v N(0_,Rp. State vector X^. is known 
to change with respect to time according to the equation 

X t = f (X^.t) + W* (2.2) 

T 

where f is some function, W_ t ^ N(0_,Q^) and E(W_ t ) = 0. 

If it is assumed that f and h in equations (2.1) and (2.2) are 
linear, then numerical procedures allow us to transform these equations 
to the form 


£k = F K ^K-l + 


= H K + Ik 


(2.3) 

(2.4) 


where K represents discrete values of t, W„ ~ N(0, Q„) and V v v 

N(0,R k ). 

After a Space Shuttle flight has taken place, the post flight reconstruc- 
tion procedure seeks to use the observed values of Z K and the equations in 
(2.3) and (2.4) to reconstruct (estimate) the parameters so that the estima- 
tion error is minimized. The estimated value of is denoted by X^. 

3. A RECURSIVE ESTIMATOR OF X 

In order to complete the development of the reconstruction process 
mentioned in section 2, we consider a system of equations 


Z = HX + e (3.1) 

where (3.1) is a system of linear equations that has been generated by 
making n observations on the single equation 

Z* - H* X + e* 


(3.2) 
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where Z* and e* are vectors or scalars. H" 1 is a jxp matrix where 
j = number of components in vector Z_* and X. is fixed. Vector X 
being fixed means that X_ does not change with time or does not change 
as the number of observations on ^increase. If Cov(t) = Z and 
det (E) r 0, then the weighted least squares estimate of •vector X is 
given by 

X = (H T e' 1 H)' 1 H T e' 1 Z (3.3) 

and the covariance of the estimate is given by 

Cov(X) = E [(X - X) (X - X) T ] = (H T E _1 H)" 1 = P . (3.4) 


The reader should be reminded that with the appropriate assump- 
tions , the maximum likelihood , the minimum variance unbiased linear as 
well as the minimum mean-squared error estimator of X are identical to 
the estimator given by equation (3.3). We also remark that Cov(X) is a 
measure of the quality of the estimates in vector X. A detailed dis- 
cussion of these estimators may be found in Elbert (1984). 

In the sequel that follows, the estimator in equation (3.3) will be 
rearranged so that X can be estimated by using a recursive process. 

The development of the recursive process is started by assuming that 
there are K observations on equation (3.2) which form the system 


— K = H K — + -K 

where Cov £ K = Z R and det ( E R ) * 0. By equation (3.3) 
= (H K T Z K 1 H K } 1 H K ^K 1 -K 


(3.5) 


(3.6) 


is the estimate of X where observations up to and including observation 
K are used. The covariance of the estimator is denoted by P v and 

P R = E[(X k - X) - X) T ] = (H k T I^ 1 h k )" 1 . 


(3.7) 
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Assume that an additional measurement Z R+1 has been made. Let 


^K+l = H K+1 — + -K+l 


(3.8) 


T T 

where Cov (£j^ + ^) = E R+1 ® -K+l^ = ^ ote ® ^£k £k+1^ 

means that the (K+l)th observation is independent of the 1st K observa- 
tions. Putting the K observations with the (K+l)th observation yields 
the system 


= 0 


ha \ 

^K+l 


\ /-K 

1 = 1-— )x + 

hr 4.1 / \ H K+1 


(3.9) 


."K+l 


By substituting into equation (3.3) we get 


&■ 


(H 


K 


H 


K+l 


-1 
K 

)l 



-rl 


£ k + x' ' H K+1' - 


(h k H K + 1> 


.-1 

‘K 


-K 


0 

S K+l' K+l 


(3.10) 


Note that we can write 


-K+l ~ + ( -K+l ’ ^K ) 


(3.11) 


Substituting results from (3.6) and (3.10) into (3.11) and using matrix 
algebra, we get 

5k+i * 5k + <p k * h k+i E ic 1 W ' 1 h k+i e k+i ( 2k+! ’ h k+i S k > • < 3 - 12 > 


By substituting into equation (3.7) and using matrix algebra, we get 


K+l 


< h kI h k+i> 


1 

'K 



-i-l 


'K+l 'H 


= (P K + h r+1 h r+1 ) 


K+l 


Substituting P R+1 into equation (3.12) yields 
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— K = -K + P K+1 H K+1 Z K+1 ( -K+l ' H K+1 ~K 3 


(3.13) 


and 


K+l 


= (P 


- 1 
K 


uT y ~ 1 U \" 

K+l K+l h K+1 j 


(3.14) 


Equations (3.13) and (3.14) provide recursive equations for estimating 
X and its corresponding covariance matrix P. 

Recall that the estimator given in equation (3.3) may be identified 
by using the symbol X^. Therefore, if X_ is fixed as it is in equation 
(3.1), the recursive equations of (3.13) and (3.14) have no statistical 
advantage over the estimator given in equation (3.3). However, if X 
changes with time where there is one observation per time interval, 
then the recursive equations are quite useful. In fact, equations (3.13) 
and (3.14) when combined with equations (2.3) and (2.4) form the 
Kalman filtering process. 


4. THE KALMAN FILTER 

It is worthwhile to mention that Brown (1983) and Gelb (1974) 
used matrix theory to express equations (3.13) and (3.14) as 


P K+1 = (I ‘ K K+1 H K+1 ) P K 

^K+l = ix + K K+1 ( -K+l ' H K+1 ^K } 
where 

K K+1 = P K H K+1 (Z K+1 + H K+1 P K H K+1 } 
and is called the gain matrix. 


(4.1) 

(4.2) 


(4.3) 


For the convenience of the reader, equations (2.3) and (2.4) are 
reprinted here. That is, 


5k+i 


= f k 5k + 2k 


(4.4) 


and 


-K+l = H K+1 -K+l + —K+l * 


(4.5) 


When equations (4.1) and (4.2) are combined with equations (4.4) and 
(4.5), two types of estimates of X R+1 are possible for each K. A 
notation by Gelb (1974) allows for the two estimators to be distinguished. 
That notation is 

X^ + ^ = the estimate of X^ + ^ using all observations up to and 
including observation K. 

= the estimate of P using all observations up to and including 
observation K. 

X^ + ^ = the estimate of using all observations up to and 

including observation (K+l). 


The values of X^^ and P^^ are obtained by using equations (4.4) 
and (4.5). The X^ +1 vector is often called the extrapolated or pre- 
dicted value of Xj£ + ^ and P^^ is called the extrapolated variance of 

X ( ' ) 

-K+l * 


The values of and P R £^ are computed by substituting X R ^ 

for X R in equation (4.2) and P K ^ for P R in equation (4.1). 


The computation summary is 
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['? = Fy Pt +) Fj + Q 


K+l 


K K r K T 


(4.7) 


V ( "^ ) — v (”).?/• cj -H V 

-K+l “ -K+l K+l ''-K+l n K+l -K+l ; 


(+) _ 
K+l 


P,/ L / = (I - K^, H v , . ) P 


(-) 

K+l “K+l' 4 K+l 


(4.8) 

(4.9) 


v — p ( ” ) u T / p p (*) uT . p n 

K K+1 ~ *K+1 H K+1 '■ h K+ 1 K+l “K+l n K+r 


(4.10) 


Equations (4.6) to (4.10) completely describe the Kalman filtering process 
when the original functions £ and h of equations (2.1) and (2.2) are 
lin.ear . 

However, when functions £ and h are nonlinear, the state vectors 

X R (K = t Q , t r ... t T ) are estimated through the application of an 

extended Kalman filtering procedure. This procedure requires a linear- 

★ 

ization of functions £ and h about some known state value X^.. The 
next paragraph provides a brief overview of the extended Kalman filter 
concept. 

Let be some known value of X and assume that AX is small. 

— J\ — 

A first degree Taylor series approximation of functions (2.1) and (2.2) 
may be 


+ AX * f(X*, K) + 


3 f 


L 8 -Jx=x k 


* • AX + W t 


and 


Z t * h(X R ) + 


3h 

TK 


L -JX=X R 


* • ax + V t 


(4.11) 


(4.12) 


* . * * 

If it is assumed that X^ is selected so that X^ = f(X^, K), then 
equations (4.11) and (4.12) become 
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AX = 


3 f 


l^JX=X K 


* “X + W K 


(4.13) 


and 



3 h 

L a -JX=X 


* AX + v K 


- -K 


(4.14) 


Equation (4.13) is called the linearized dynamics equation and (4.14) 
is the linearized measurement equation. Note that equations (4.13) and 
(4,14) may be transformed to equations that are equivalent to equations 
(2.3) and (2.4). Hence for each discrete time K, aX k can be estimated 
and error covariance matrices can be determined. The state vector 
estimate at time K is then given by 


5k +) * *K 


A X^ 

-K 


where aX r is computed by substituting from equations (4.13) and 

(4.14) into equation (4.8). Vector X K+1 is computed by letting 
* 

-K+l = -K anc * re P eat i n ? the above procedure. 

The extended Kalman filter has performed well in a large class of 
applications. However, there are occasions when divergence occurs in 
the state vector estimates. Divergence occurs when the computed 
entries of the error covariance matrix P K become small as compared to 
the actual error in the estimate of the state vector. The cause of this 
divergence is not due to a defect in the filtering procedure, but may 
be caused by the linear approximation procedure, numerical rounding 
error, or an inadequate model of the system being studied. Additional 
possible causes of divergence are mentioned by Gelb (1974). 



A publication by Varhaegen and Van Dooren (1986) is representa- 
tive of the theoretical and experimental analyses that are currently 
being done on the divergence problem. The Rogers (1987) computer 
code has implemented the U-D factorized algorithm as a means of con- 
trolling numerical roundoff error that may lead to divergence in the 
state vector. Checks for other sources that may generate divergence 
are conducted through computer evaluations. 

5. CONCLUSION 

At this point we have reviewed the general setup for a regression 
problem where the parameters to be estimated are fixed. This lead to 
recursive equations (4.1), (4.2), and (4.3) which allowed for the 
development of an estimation procedure for a time varying parameter. 
When functions £ and h are linear, it has been clearly stated that 
equations (2.3) and (2.4) are the essential ingredients for state 
variable reconstruction. If f and h are nonlinear, can be estimated 
by linearing h and f about some known vector and treating the 
linearized equations as if they are equations (2.3) and (2.4). 

The problem of applying the extended Kalman filtering procedure 
to the static test firing data remains. The basic approach for the 
application is identical to the presentation given in Section 4. There- 
fore, the computer codes that have been prepared by Rogers (1987), 
which are currently operational on the MSFC computer system, may be 
modified so that static test data may be analyzed. The static test 
setup differs from the actual flight data in that many flight associated 
modules will become inactive. Therefore, after adjustments are made 
for the analysis of static test data, it will be necessary to evaluate the 
performance of the computer code for the divergence of parameter 


estimates. 
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